A MAPPING FROM THE UNITARY TO DOUBLY 
STOCHASTIC MATRICES AND SYMBOLS ON A 

FINITE SET 

ALEXANDER V. KARABEGOV 

Abstract. We prove that the mapping from the unitary to the 
doubly stochastic matrices that maps a unitary matrix (v,ki) to 
the doubly stochastic matrix (|wfc;| 2 ) is a submersion for almost 
all unitary matrices. The proof uses the framework of operator 
symbols on a finite set. 



1. Introduction 

Given a unitary matrix (%), the matrix (|mh| 2 ) is doubly stochastic. 
The mapping 

fJ> : (Uki) (Kz| 2 ) 
from the set U of unitary (n x n)-matrices to the set P of doubly 
stochastic (n x n)-matrices was studied in a number of papers in math- 
ematical physics (see [5] and references therein). The set U is a real 
n 2 -dimensional manifold and P is a convex (n — l) 2 -dimensional poly- 
tope, the Birkhoff polytope. It is known that the image fi(U) is a 
proper subset of the Birkhoff polytope P. Two unitary matrices u and 
u are called equivalent if there exist unitary diagonal matrices x and A 
such that u = xu\. The mapping \x maps equivalent unitary matrices 
to the same doubly stochastic matrix. The set of equivalence classes of 
unitary matrices with nonzero entries is an (n — l) 2 -dimensional mani- 
fold. These dimensional considerations lead to the natural question of 
whether \i is a submersion. It was deemed that \i has to be a submer- 
sion almost everywhere on U (see, e.g., [T]). In this paper we give a 
detailed proof of this result which was announced in [3]. 

To each unitary (n x n)-matrix u with nonzero entries we relate 
an operator I u on an n 2 -dimensional vector space V and show that 
the dimension of the kernel of the tangent mapping equals the 
multiplicity of the eigenvalue 1 in the spectrum of I u . We introduce 
the notion of a symmetry group G of a unitary matrix u. The group 
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G has a unitary representation in the space V which commutes with 
the operator I u . When the matrix u has a large symmetry group, 
it is sometimes possible to calculate the spectral decomposition of the 
operator I u . This allows us to show that the mapping /j, is a submersion 
for almost all elements of U. 

2. Symbols on a finite set 

In this section we introduce a construction of two types of symbols 
of operators on a finite dimensional vector space. These symbols are 
analogues of pq- and gp-symbols in quantum mechanics (see [2]). 

We will model finite dimensional vector spaces as the spaces of func- 
tions on finite sets. Given a finite set J, denote by F(J) the space of 
all complex-valued functions on J. The standard Hermitian product of 
functions (p — (pj,ip — ipj G F(J) is given by the formula 

(<f,ip) = J^^r 

Now let K and L be two n-element index sets and u = (uki),k G 
K, I G L, be a unitary matrix. It determines a unitary isomorphism 
U : F(L) —>■ F(K) of the vector spaces F(K) and F(L) endowed with 
the standard Hermitian products: for -0 — ipi G F(L), 

leL 

Denote M := K x L. We define two mappings 

C U ,D U : F(M)^ End F(K) 

such that for / = /jy G F(M) the matrices (xkk>) and (ykk 1 ) of the op- 
erators C u f and D u f in F(K), respectively, are given by the following 
formulas: 

(1) %kk' = }u k if k iu k 'i and y kk > = ^2u k if k 'iU k n. 

leL leL 

The mappings / i— ► C u f and / i— > D u f can be thought of as two 
different symbol-to-operator mappings. 

Given functions a = a k G F(K) and b = bi G F(L), denote by d and 
b the multiplication operators by these functions in F(K) and F(L), 
respectively. Then ab = a k b\ G F(M), UbU* G End F(K), and it can 
be checked that 

(2) C u (ab) = a (ubU*) and D u (ab) = (ui>U*) a. 
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Using the analogy with quantum mechanics, we may think of F(K) 
and F(L) as of the "coordinate" and "momentum" representations of 
an abstract Hermitian vector space, respectively, with the isomorphism 
U playing the role of the Fourier transform. The mapping C u is ob- 
tained by the standard ordering of the "coordinate" operators a and 
"momentum" operators UbW, while D u is obtained via the inverse 
ordering. 

Given a function / e F(M) and an operator A e End F(K), we will 
say that the function / is a C-symbol of the operator A if A = C u f 
and that / is a D-symbol of A if A = D u f . Thus, C- and D-symbols 
correspond to qp- and pg-symbols, respectively. 

Now we will make a critical assumption that all elements u^i of the 
matrix u are nonzero and introduce a Hermitian product (•,•)« on the 
space F(M) by the formula 



The space End F(K) of operators in F(K) carries the Hilbert-Schmidt 
Hermitian product (•, -)hs- If ( x kk') an d (ykk') are the matrices of 
operators X and Y in F(K), respectively, then 



The following properties of C- and D-symbols are analogous to the 
corresponding properties of pq- and gp-symbols of Hilbert-Schmidt op- 
erators in quantum mechanics. 

Proposition 2.1. The mappings C u , D u : F(M) — > End F(K) are uni- 
tary isomorphisms of the Hermitian vector spaces (F(M), (•, -) u ) and 



Proof. Given two functions f,g e F(M), let (x^/) and {ykk 1 ) be the 
matrices of the operators C u f and C u g, respectively, so that 




keK,ieL 



(X,Y) HS = trXY* = ^ x kk'Vkk'- 



k,k'eK 



(EndV,(-,-) HS ). 




Then, using the unitarity of the matrix ti, we obtain that 
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where 5w is the Kronecker symbol. Thus C u : F(M) — > End F(K) is 
an isometric mapping. Since the dimensions of F(M) and End F(K) 
are both equal to n 2 , the mapping C u is a unitary isomorphism. A sim- 
ilar calculation shows that D u : F(M) — > End F(K) is also a unitary 
isomorphism. □ 

It follows from Proposition 12.11 that the mapping I u := C~ l D u that 
maps the D-symbols to the corresponding C-symbols is a unitary op- 
erator in (F(M), (•, •}„). Using the analogy with quantum mechanics, 
we will call the operator I u the Berezin transform. It can be checked 
that the operator I u is given by the following formula: for a function 

f = f kl eF(M), 

I J f\ u kl'Uk'l f i ,2 

\J-uJ)kl— / J Jk'l'\Uk'l'\ ■ 

z — ' UklUh'l' 

k'eKd'eL M K 1 

Proposition 2.2. Given a junction f G F(M), the operators C u f and 
D u f are Hermitian conjugate. 

Proof. Let {x^y) be the matrix of the operator C u f in F(K). Using 
formulas ([1]) we see that the Hermitian conjugate matrix (x^k) is given 
by the equation 

Xk'k = 'y^ J Uk'lfk'lUkl 

and therefore coincides with the matrix of the operator D u f. □ 

This proposition allows to give a characterization of the C- and D- 
symbols of the skew- Hermitian operators in F(K) which will be used 
in Section 

Corollary 2.1. Given functions f,g £ F(M), the operator C u f is 
skew- Hermitian if and only if f = —I u f an d the operator D u g is skew- 
Hermitian if and only if I u g = —g. If f = I u g, then the operator 
C u f = D u g is skew- Hermitian if and only if f = —g. 

Proof. Proposition 12.21 states that 

(cjy = dJ. 

Now the operator C u f is skew-Hermitian, i.e., C u f = —(C u f)*, if and 
only if C u f = —D u f or, equivalently, / = —I u f. The rest of the 
corollary is proved similarly. □ 

Assume that functions /, g G F(M) are such that / = I u g or, equiv- 
alently, C u f = D u g. According to Proposition 12.21 

Dj=(C u fr = (D u g)* = C u g, 
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whence I u f = g. Let g G F(M) be an eigenfunction of the Berezin 
transform I u with the eigenvalue 9. Since I u is unitary, we have \9\ = 1. 
Equation I u g = 9g implies that I u {9g) = g, whence I u (g) = 9g. Thus 
we have proved the following 

Lemma 2.2. The eigenspaces of the Berezin transform I u are invari- 
ant under complex conjugation. In particular, for any eigenvalue 9 of 
multiplicity k the 9-eigenspace treated as a real vector space is a di- 
rect sum of the k-dimensional real vector spaces comprised of real and 
purely imaginary eigen] 'unctions, respectively. 

Denote by E the subspace of F(M) of the functions / that can be 
represented as a sum of the form fki = at + h. The dimension of E 
is 2n — 1. Formulas (|2J) imply that for functions a = and b = bi 
considered as elements of F(M), 

C u (a) = D u (a) and C u (b) = D u (b), 

which leads to the following statement. 

Lemma 2.3. For any function f G E we have I u f = f . Thus the 
multiplicity of 1 in the spectrum of the Berezin transform I u is at least 
2n - 1. 

3. Symmetries of a unitary matrix 

Assume that, as in Section^ K and L are two n-element sets, and 
u — (uki), k G K, I G L, is a unitary matrix which defines a unitary 
isomorphism 

U:F{L)^F{K), 

where F(K) and F(L) are endowed with the standard Hermitian prod- 
ucts. Assume further that G is a group that acts on the left on the 
sets K and L and has faithful unitary representations S and T in the 
spaces F(K) and F(L), respectively, such that, given g G G,ip = ifk G 
F(K),i> = i> l eF(L), 

{S g (f) k = a k (g)(f g -i k and (T g ip) l = bi(g)if)g-i h 

where a k (g) G F(K) and bi(g) G F(L) are two unitary functions (i.e., 
^(fiOl — \bi(g)\ = 1)- Finally we assume that the isomorphism U 
intertwines the representations S and T: for g G G, 

(3) S g U=UT g . 

If these conditions are satisfied, we say that G is a symmetry group of 
the matrix u. Condition ([3]) can be expressed in terms of the matrix u 
as follows: 



(4) 



a k{g) {ug-ikg-n) bi(g) = u M - 
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Denote by R the following representation of the group G in F(M): for 
geG and f = f kl eF(M), 

( R gf)ki = fg-^kg-H- 
Using formula (BJ and the fact that the functions a k (g) and bi(g) are 
unitary we obtain the following 

Lemma 3.1. The representation R is a unitary representation of the 
group G in the Hermitian vector space (F(M), (-, •)„). 

Denote by Q the representation of the group G in the space End F(K) 
by conjugation via the operators of the representation S: for g G G 
and X G End F(K), 

Qg(X) = SgXS g 1 = SgXSg. 

The representation Q is a unitary representation of the group G in the 
Hermitian vector space (End F(K), (-, -)hs)- 

The following equivariance property of the symbol-to-operator map- 
pings C u and D u is analogous to the equivariance of the pq- and qp- 
symbol mappings with respect to the Heisenberg group. 

Proposition 3.1. The mappings C U ,D U : F(M) — > EndF(K) inter- 
twine the representations Q and R: for g G G, 

C u (R g f) = Sg(C u (f))S; andD u (R g f) = S g (D u (f))S* g . 

Proof. Given g G G, f G F(M) and ip G F(K), we get with the use of 
formula (jlj) the following: 

(SC u (f)ip)k = a k(g)u g -iki> fg-iki'Uk>i"Pk' = 

k'£K,l'eL 

a k{g)Ug-Xk9-Hf9-ikg-HU g -r kg -i l <p g -i~ k = 

k£K,ieL 

k£K,l£L 

U ^ R 9 f)klHg)Ug-, k g- n ^g-l k = 

k£K,leL 

u ki{R a f)ki (a k {g)u g -i kg -iM9)) a kig)Va{k) = 

k£K,ieL 

u ki(Rgf)kiu kl a k (g)p g -i k = (C u (Rgf) S<p) k . 

k£K,leL 

The proof that the mapping D u intertwines the representations Q and 
R is similar. □ 
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Since both C u and D u intertwine the representations Q and R of the 
group G, we obtain the following 

Corollary 3.2. The Berezin transform I u = C~ 1 D U commutes with 
the representation R: for g £ G, 

4. Examples of unitary matrices with a large group of 

symmetries 

In this section we will consider two examples of unitary matrices 
with large symmetry groups. For each of these matrices we will obtain 
the spectral decomposition of the corresponding Berezin transform and 
in particular determine the multiplicity of the eigenvalue 1 in their 
spectra. In Section [5] we will give a geometric interpretation of this 
multiplicity. 

Example 1. Fix a natural number n and denote by e the character 
of the cyclic group Z/nZ given by the formula 

f 2mk \ 



e{k) = exp 



Here we denote by k (somewhat loosely) an integer modulo n. The 
Fourier transform T on the cyclic group Z/nZ is given by the unitary 
matrix u = (uki) with 

u k i = —j=e(kl), 
V n 

where k,l £ Z/nZ (here K = L = Z/nZ). The corresponding Berezin 
transform is given by the following formula: 

(5) = ^ Yl e(-(k-k')(l-l'))f kn/ . 

k',l'eZ/nZ 

For r £ Z/nZ denote by W r the unitary operator of multiplication by 
the function e(rk) and by Z r the shift operator by the element r such 
that for ifk £ F(Z/nZ) 

{W r Lp) k = e(rk)ip k and (Z r ip) k = <p k +r- 
A simple calculation shows that 

(6) T*W r T = Z_ r and T*Z r T = W r . 
Formulas (j2J) and (jSJ) imply that 

C u (e(rk + si)) = W r TW s T* = W r Z s . 

Similarly, 

D u (s(rk + si)) = Z s W r . 
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From the Weyl commutation relations 
(7) Z s W r = e(rs)W r Z ! 



s 



it follows immediately that 

I u (e(rk + si)) = e(rs)e(rk + si). 

The functions e(rk + si) with r, s G Z/nZ form a complete system 
of eigenvectors of the operator I u with the corresponding eigenvalues 
e(rs). Observe that e(rs) = 1 if and only if rs = in Z/nZ, which 
implies the following 

Lemma 4.1. The multiplicity of the eigenvalue 1 of the Berezin trans- 
form corresponding to the matrix of the Fourier transform on the group 
Z/nZ is equal to the number of pairs (r,s) of integers such that < 
r,s < n and rs = (mod n) . In particular, this multiplicity equals 
2n — 1 if and only if n is a prime. 

Remark. The spectral decomposition of the Berezin transform I u is easy 
to obtain because the matrix u of the Fourier transform on the group Z/nZ 
has a large symmetry group. Let T denote the multiplicative group of 
complex numbers of modulus 1. The generalized Heisenberg group G := 
(Z/nZ x Z/nZ) X T with the product 

(si, t\, 0i) • (s 2 ,t 2 ,6 2 ) = (si + S2, t\ + t 2 , 9i0 2 e(hs 2 )), 

where Si,ti £ Z/nZ and 6i £ T, is a symmetry group of the matrix u. The 
corresponding representations S and T are given by the formulas 



S(s,t,e) = 0W s Z t and T {sm = F*S (tm F = 6Z_ s Wt = 8e(-si)W t Z_ 



Given functions (fk and V>z on Z/nZ and an element (s,t,6) G G, we have 



Since the Berezin transform © commutes with the representation R it is 
shift-invariant and thus the functions e(rk + si) with r,s G Z/nZ are its 
eigenfunctions. 

Example 2. Let n be a fixed natural number. Set K = {1, 2, • • • , n} 
and denote by ©„ the group of permutations of K (the symmetric 
group). Fix a constant 9 such that |#| = 1 and 6 ^ ±1. A simple check 
shows that the matrix w = (um), k,l & K, with 

^ 

(8) u H = 5 H H 
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has nonzero elements and is unitary. Denote by U G End F(K) the 
corresponding unitary operator (we assume that F(K) carries the stan- 
dard Hermitian structure). Let S be the unitary representation of the 
symmetric group & n in the space F(K) by permutations: for a G & n 
and if = ipk G F(K), 

Since the operators U and S a commute, the group & n is a symmetry 
group of the matrix u. According to Lemma 13.11 the representation R 
of the group & n in the Hermitian vector space (F(M), (■,■)«) given by 
the formula 

is unitary. Here a G & n and / = /jy G F(M). 

The representation R has the following four isotypic components. 
Denote by V\ the two-dimensional space of functions from F(M) of 
the form f ki = a + b8ki, where a and b are constants. The group & n 
acts trivially on V\. The space V2 C F(M) of functions of the form 
fki = a k + bi + c k 5 M where 

n n n 

y]°< = & » = 5Z q = °' 

j=l i=l j=l 

is the isotypic component of multiplicity 3 of the representation R 
corresponding to the irreducible representation of G n of dimension n — 1 
parametrized by the partition (n — 1,1) according to [3|. The space 
V3 C F(M) of functions f k i such that f k i = —fik f° r all k,l <E K and 

n 

(9) X> = ° 

is irreducible. It is parametrized by the partition (n — 2, 1, 1) and has 
the dimension 

, Tr n 2 -3n + 2 (n-l)(n-2) 

dim V3 = = . 

2 2 

Finally, the space V4 C F(M) of functions such that /^z = fi k , f kk = 
for all k,l E K, and satisfying condition (jHJ) is also irreducible. It is 
parametrized by the partition (n — 2, 2) and has the dimension 

n 2 — 3n n(n — 3) 

dim Va = = . 

2 2 

The space E C F(M) of functions of the form f k i = a k + bi is & n - 
invariant. Therefore, its orthogonal complement E 1 - with respect to 
the Hermitian product (•,•)« is also (5 n -invariant. The restriction of 
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the representation R onto E x is multiplicity free. It is the following 
orthogonal sum of irreducible subspaces: 

(io) E 1 = (Vi n E 1 ) © (v 2 n e l ) © v 3 © v A . 

Using the facts that the Berezin transform I u commutes with the repre- 
sentation R and acts trivially on the space E, we will find the spectral 
decomposition of the operator I u . It is scalar on the direct summands 
in (TIP]). 

Given a function f = fki & V3, we will calculate the matrices (xkw) 
and (ykk 1 ) of the operators C u f and D u f, respectively. Observe that 
fkk = for all k E K. Using formulas ([!]) and Qj, we get 



Xkk' 



£ 

lei 



leL 
1 



0-1 



it; 



7? 



A7 



Sklfkl^k'l H /fcA'Z + 5jw/fcJ 

n n 



0- 



0-1 



7? 



77 



/i 



A7 



fkk^kk' + 







7? 



■Jfefe' + Jkk- 



n 



?? 



-fi 



kk J 



A similar calculation shows that 







ykk' 



fk'k 







/i 



A' A' 



/fefc' 



n 



£fcfc' • 



n n 

Thus D u f = 6 C u f, whence I u f = Of. 

Now assume that / = fw G V4. Denote by and (y^/ 

matrices of the operators C u / and respectively. Then similar 

calculations show that 



the 



0-1 



Xkk' 



fkk' and 







A; A' 



-Oxkk' 



n n 

whence D u / = —0C u f and therefore I u / = —0/. 

Tedious but straightforward calculations render the eigenvalues of 
the operator I u on the eigenspaces V\ fl E 1 - and V2 PI -E -1 . The spec- 
tral decomposition of the Berezin transform I u is summarized in the 
following table. 



Eigenspace 


E 


^nE 1 


v 2 nE L 


V 3 


u 4 


Dimension 


2n- 1 


1 


n — 1 


n 2 -3n+2 
2 


n 2 — 3n 
2 


Eigenvalue 


1 


ne+n-l 
° 9+n-l 


0+n-l 
0+n-l 





-0 



Observe that the condition 7^ ±1 implies that the multiplicity of 1 
in the spectrum of the Berezin transform I u is equal to 2n — 1. 
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5. A MAPPING FROM THE UNITARY TO THE DOUBLY STOCHASTIC 



Let K and L be two n-element index sets. Denote by U the set of 
all unitary matrices (uki), k G K,l G L. It is a principal homogeneous 
space of the group of unitary matrices (ukk'),k,k' G K, and thus is 
a real n 2 -dimensional compact manifold. Denote by Uq the set of all 
unitary matrices from U with nonzero entries. It is an open subset of 



Let B be the affine space of real matrices (pu), k G K,l G L, satis- 
fying the equations 



Its dimension is (n — l) 2 . The set P of matrices from B with nonnega- 
tive elements is an (n — l) 2 -dimensional convex polytope, the Birkhoff 
polytope (see |6j). Its elements are called doubly stochastic matrices. 
The Birkhoff polytope is the convex hull of the permutation matrices 
(^cr(fc)z), where a : K — > L is a bijection. Its interior consists of the 
doubly stochastic matrices with nonzero entries. 

Two unitary matrices u = {uki) and u = (ukl) from U are called 
equivalent if there exist unitary diagonal matrices x = {xk&kk 1 )-, k, k' G 
K, and A = (XiSw),l,l' G L, such that u = >tu\, i.e., Uki = ^kUkAi- 
We will write u ~ u for equivalent matrices u and u. 

Denote by // : U — > B the mapping that maps a unitary matrix (uki) 
to the doubly stochastic matrix (|-Ufcz| 2 ). It is known that the image 
n(U) is a proper subset of the Birkhoff polytope P (see [5]). 

The mapping /x maps equivalent matrices to the same doubly sto- 
chastic matrix. Thus it factors through the singular qoutient space 
Uf ~. The quotient space Uq/ ~ is an (n — l) 2 -dimensional manifold. 
These observations lead to the natural question of whether the map- 
ping fj, is a submersion. It is obvious that \i is not a submersion on the 
complement U\Uq since it is mapped to the boundary of the Birkhoff 
polytope. In the rest of this paper we will show that for each n there 
exists a unitary (n x n)-matrix u such that \i is a submersion at u. 
Since /i is an algebraic mapping of real affine manifolds, this will imply 
that fj, is a submersion for almost all elements of U. 

Let it = (um) be an arbitrary unitary matrix in U and U : — > 
F(K) be the corresponding unitary isomorphism. We want to give 
different descriptions of the tangent space TJJ . Assume that u(t) = 
Uki(t) G U is a family of unitary matrices smoothly depending on a 
small real parameter t and such that it(0) = w. Then u'(0) = (u' kl (0)) 
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is a tangent vector in T U U. Here we use that U is a real submanifold 
in the affine space of complex- valued matrices (xki) with k G K and 
/ G L. This way one can obtain all tangent vectors in T U U. Each 
matrix u(t) determines a unitary isomorphism U(t) : F(L) F(K). 
Since U(t)(U(t))* = 1 and W(0) = W, we get that 

U'(0)U* +U(U'(0))* = 0. 

The operator 

(11) x = u'(o)u* = -u(u'(o)y 

in F(K) corresponding to the tangent vector tt'(0) is skew-Hermitian. 
This correspondence identifies the tangent space TJJ with the space 
of skew-Hermitian operators in F(K). 

Now assume that u G Uq. Corollary 12 . 1 1 allows to identify the space of 
skew-Hermitian operators in F(K) (and thus the tangent space T U U$) 
with the space of C-symbols / G F(M) such that / = —I u f and also 
with the space of .D-symbols g G F(M) such that I u g = —g. 

Formula fill I) implies that the matrix (xkk'), k, k! G K, of the operator 
X satisfies the following equations: 

(12) X k¥ = ^u' kl (Q)u kn = -^2u kl u' kn (0). 

leL leL 

Assume that / = f k i and g = g k i are the C- and D-symbol of the 
operator X, respectively, that is, 

X = C u f = D u g. 

It follows from formulas (|T]) and ffl2l) that 

Ukl u H 
In particular, / = —g. Formulas (fTTj) and (fl3|) mean that the C- and 
D-symbols are related to the perturbations of the matrix elements of u 
while the operators corresponding to these symbols are related to the 
perturbations of the isomorphism U. 

The mapping \i maps the matrix u{t) = (u k i(t)) to the doubly sto- 
chastic matrix p(t) = (p k i(t)) such that 

Pkl{t) = \Ukl{t)\ 2 = U k i(t)Ukl{t). 

We set p = p(0). The tangent mapping maps the tangent vec- 
tor u'(0) G T u Uq to the tangent vector p'(0) G T p P . It follows from 
formulas (ITS]) that 

P«(°) = u 'ki(fyuki + u k iu' H (0) = (fki ~ 9m) \uh\ 2 - 
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Thus the tangent vector tt'(O) G T u lIo lies in the kernel of the tangent 
mapping \i^ u if and only if the C- and D-symbols of the corresponding 
skew-Hermitian operator X coincide, / = g, or, equivalently, I u g = g, 
that is, g is an eigenvector of the operator I u with the eigenvalue 1. 
Since the operator X is skew-Hermitian, according to Corollary 12. 1\ 
f = —g and therefore the condition that f = g implies that g = —g, 
i.e., that g is purely imaginary. We conclude that the kernel of the tan- 
gent mapping can be identified with the space of purely imaginary 
eigenfunctions of the Berezin transform I u with the eigenvalue 1. Now, 
Lemma [2.21 implies the following 

Theorem 5.1. Given a unitary matrix u G Uq, the dimension of the 
kernel of the tangent mapping fi* u is equal to the multiplicity of 1 in 
the spectrum of the Berezin transform I u . 

Theorem 15.11 and Example 2 from Section H] show that the mapping 
fj, is a submersion at the matrix u = {uu) given by (jSJ). Thus, \i is a 
submersion almost everywhere on U. 
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